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Abstract. The vacuum solutions arising from a spontaneous breaking of Lorentz symmetry due to the acquisition of a 
vacuum expectation value by a vector field are derived. These include the purely radial Lorentz symmetry breaking (LSB), 
radial/temporal LSB and axial/temporal LSB scenarios. It is found that the purely radial LSB case gives rise to new black hole 
solutions. Whenever possible, Parametrized Post-Newtonian (PPN) parameters are computed and compared to observational 
bounds, in order to constrain the Lorentz symmetry breaking scale. 



1, INTRODUCTION 

Lorentz invariance is clearly one of the most fundamental symmetries of Nature. It is both theoretically sound and 
experimentally well tested [1, 2], thus playing a leading role in most theories of gravity. Therefore, it is only natural 
that little attention has been paid to the consequences of explicitly breaking this symmetry. 

A more flexible approach to this question admits a spontaneous breaking of this symmetry, instead of an explicit 
one, analogously to the Higgs mechanism in the Standard Model of particle physics [3, 4, 5, 6]. This can arise if a 
vector field ruled by a potential exhibiting a minimum rolls to its vacuum expectation value (vev) - this vector field, 
usually referred to as "bumblebee" vector, thus acquires a specific four-dimensional orientation. 

From a theoretical standpoint, a spontaneous Lorentz symmetry breaking (LSB) is possible, for instance, in 
string/M-theory, arising from non-trivial solutions in string field theory [3, 4] and in noncommutative field theories 
[7, 8]. A spacetime variation of fundamental coupling constants could also lead to a spontaneous LSB [9]. Experimen- 
tally, the violation of Lorentz invariance could be tested in ultra-high energy cosmic rays [10]. 

The consequences of the "bumblebee" vector scenario were studied in Ref. [11]; in there, three relevant cases were 
taken into account: the bumblebee field acquiring a purely radial vev, a mixed radial and temporal vev and a mixed 
axial and temporal vev. The results were analyzed in terms of the PPN parameters, when possible, prompting for 
comparison with current and future experimental bounds and effects, for instance, from string theory in a low-energy 
scenario. These bounds may arise from the observations of the Bepi-Colombo [12] and LATOR [13] missions (see 
also Ref. [14]) for a discussion on future gravitational experiments). 

The action of the bumblebee model is written as 



(1) 



where k = %kG, B^v — d^By — dyB^, <^ is a coupling constant and sets the bumblebee's vev, since the potential 
V driving Lorentz and/or CPT violation is supposed to have a minimum at B^B^ ±b^ —0, i.e. V'{b^b^) = 0. The 
particular form of this potential is iiTelevant, since one assumes that the bumblebee field is frozen at its vev. The scale 
of b^ should be obtained from string theory or from a low-energy extension to the Standard Model. Hence, one expects 
bj^ to be of order of Mpi, the Planck mass, or Mew, the electroweak breaking scale. 



2. PURELY RADIAL LSB 



In this section, a method to obtain the exact solution for the purely radial LSB is developed. A static, spherically 
synunetric spacetime, with a Birkhoff metric gfiv = diag{—e^^,e^f,t^,r^sirP'0) is considered. It can be easily seen 
that the Killing vectors of the metric are conserved, showing that radial symmetry is still valid; this enables the 
construction of a covariantly conserved current associated with the vector bumblebee field [11]. 

The affine connection derived from the metric g^y allows for the computation of b^i, given the non-trivial covariant 
derivative with respect to the radial coordinate, and taking fo^ = (0, /^(r) , 0, 0) . Hence, from D^bv = d^bv— Tv/j^a = 0, 

it follows that b{r) = £,~^l'^b(,eP , where the factor is introduced for later convenience. As expected, b^ = b^b^ = 
bl^~^ is constant. 



The (spatial) action can be thus written as 



/ 



R 



- + (/'-)V(r)/?. 



(2) 



The determinant of the metric is given by = P'eP^^; the scalar curvature and the relevant non-vanishing Ricci 
tensor component are given by = 2r"^ [l -|- (2rp' — 1 )e"^''] and Rrr = 2r^^p', where the prime stands for derivative 
with respect to r and we have integrated over the angular dependence. Also, ^ {b'')-^Rrr = bQ2r~^p'e~^P, where ^ is 
the contravariant radial component of b^. By introducing the field redefinition I* = (l — e~^f) r~^, the action may be 
rewritten as [15] 



Variation with respect to produces the equation of motion 



(3 + fc2)^+(l + J)y^' 



(3) 



{3+blYV+{+l + -^)r^' =Q , 



(4) 



which admits the solution »P(r) = ^or-^+^, with L = 3 - (3 -h fco)/(l + ^1/2) ~ bl/2, and hence L ~ by 2. We thus 
obtain grr = e^P = (l — I'o''"^^^) ^ • Comparison with the usual Schwarzschild metric yields 



2Gimr 

500 = -H 



(5) 



where Gl has dimensions [Gl] = L}~^ (in natural units, where c = h= \); one can define Gl = Gr^^, where ro is 
an arbitrary distance. The L ^ Umit yields Gl^ G and the usual geometrical mass, M = Gm, with dimensions of 
lenght. From now on we express all results in terms of M. The event horizon condition is given by 



goo ■■ 



2M rL 

-l-h : 
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thus rs = {2MrQ ^)^/^ ^. The norm-square of the Riemann tensor, I = R^ivpi^^^^^^ ^ is found to be 



7 = 48 
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(6) 



(7) 



where /q = 48M^r~^ is the usual scalar invariant in the limit L — > 0. Since I{r = r^) is finite, the singularity at r = 

is removable; accordingly, the singularity at r = is intrinsic, as the scalar invariant diverges there. One concludes 
that an axial LSB gravity model admits new black hole solutions with a singularity well protected within a horizon of 
radius r^. The associated Hawking temperature is 



n 1 



(2Mro^)- 



)ro-(2Mro^)-^ro 



where 7b = h/SnksM is the usual Hawking temperature, recovered in the Umit L — > 0. 



(8) 



Since the obtained metric cannot be expanded in powers ofU = M/r, a PPN expansion is not feasible. However, a 
comparison with results for deviations from Newtonian gravity [16], usually stated in terms of a Yukawa potential of 
the form 

y,(.) = ^(l + a.-A) , (9) 
yields Gir^ = Gy (l + ae~''^^), which, to first order around r = ro, reads 

Gz.^(^l+L^^ =Gr(l + a-a^) , (10) 

so that one identifies A = ro and a = —L (with Gy(l — L) = Gifo = G). Planetary tests to Kepler's law in Venus 
indicate that A = ro = 0.723 Af/ and L = |a| < 2 x 10"^. 



3. RADIAL/TEMPORAL LSB 

We consider now the mixed radial and temporal Lorentz symmetry breaking. As before, it is assumed that the 
bumblebee field B^j has relaxed to its vacuum expectation value. Provided that one takes temporal variations to be 
of the order of the age of the Universe H^^ , where //q is the Hubble constant, a Birkhoff static, radially symmetric 
metric g^y = diag{—e^'^,e^P,r^,r^sin^{9)) may still be used. The physical gauge choice of a vanishing covariant 
derivative of the field 5^ yields br{r) = ^^^I^ArC^ and, similarly, bQ{r) = B,~^I^AQe^ , with Aq and A^ dimensionless 
constants. As before, b^ = b^b^i = (A^ — Aq)(^"' is constant. 

In the present case, the symmetry = — p does not hold, as now both a radial and a temporal component for the 
vector field vev are present; for this reason, the previous spatial action formalism depicted in Eq. (3) cannot be used. 
Instead, the full Einstein equations must be dealt with. 



1 



{b^'fRaaguV - b^.h'R^y - bybf'Ry^ 



(11) 



Since the bumblebee field has relaxed to its vev and therefore both the field strength and the potential term vanish, 
the additional equation of motion for the vector field is trivial. The metric Ansatz and the expressions for b^ then yield 



1 

Goo =2 



3AlRoo-Ay^'^-P^Rr 



Ay^P-^^R^-2,A%. 



(12) 



Writing Goo = go{r)e^^^-P\ Grr = gr{r), Roo = /o(Oe2(^-P) and Rrr = Mr), where 



r r 



-l+e^P 2p' Id)' 
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the Einstein equations read 

go{r) = \ [3A2/o(r) -Alfr{r)\ , gr{r) = i [AlMr) - SA^M] • (13) 

This is the set of coupled second order differential equations which must be solved, with boundary conditions given 
by (j)(oo) = p(oo) = (j)'(oo) = p'(oo) = 0. 

The spontaneous LSB is clearly exhibited; as can be noticed from gQ+ g^ = fo+ fr, one has (1 — 2Aq)/o = 
— (1 + 2A^)fr; in the unperturbed case Ao = Ar = 0, /o = — fr, and the Schwarzschild solution (j) = —p is recovered 
from go + gr = 0. This symmetry does not hold in the perturbed case, which produces /o w — (1 + ^A^ + lA^jfr. 

An expansion of the metric in terms of (j) = (j)o + d(l) and p = — <^ — 5p allows for the solving of Eqs. (13), where 
is given by the usual Szcharzschild metric, ^ = jln (1 — 2M/r), and 5p, 50 are assumed to be small perturbations. 
After some algebra, the solution is found to be [1 1] 



where A = Aq, B = A^, is an integration constant and 



l + a^— — - 



Kr-" , (14) 



^^-C^-C^WK±M±^^0, (15) 



with 



Ci =A + 3fi+Afi + 9fi^~A + 3B , C2 = 2 + B- 3A + 16AB~ 2 , C3=2 + B-7A~2 . (16) 

One can linearize the exponent a around Ci <C 1, yielding a ~ (2 — 7A +fi)/(2 + A + 35), so that a ~ 1. 
After solving the coupled differential Eqs. (13), the non-trivial components of the metric now read 



g-2(.fc+5p) 



r 

-(2+a(A+B))/fr-« g-^^^rr" 



with the definition /T^ = [1 + a(A -|-fi)/2]A' ~ [1 + (A -|-fi)/2]A'. Following the algebra of a Lorentz transformation to a 
isotropic coordinate system, on which all spatial metric components are equal, and then to a quasi-cartesian referential, 
the resulting metric is [1 1] 

/ K + Kr\ , / K + 2Kr\ 

%=^« = -l+2f/- ( , 77^,^, = 1 + h - j f/ , (18) 

and the PPN parameters may be directly read, yielding j3 ~ 1 — {K + Kr)/M and 7~ 1 — {K + 2Kr)/M. Inverting 
this relation gives K/M ~ 1 — 2/3 + 7, Kr/M ~ j3 — 7. Hence, a temporal/radial LSB manifests itself linearly on the 
PPN parameters /3 and 7. A caveat of these results is the clear dependence of the obtained PPN parameters on the 
free- valued integration constants K and Kr, instead of the physical parameters associated with the breaking of Lorentz 
invariance. This reflects the linearization procedure followed in order to obtain the radially symmetric Birkhoff metric 
solution to the Einstein equations. 

The bounds derived from the Nordvedt effect, |j3 — 1| < 6 x 10""^ [17] and the Cassini-Huygens experiment, 
7= 1 -I- (2.1 ±2.3) X 10-5 [18], can be used to obtain \K + Kr\ < 0.9 m md K + 2Kr = (-3.1 ±3.4) x lO'^ m. 
Since, by definition 



Kr 



(A+B)' 



K , (19) 



with a ~ 1, A,B <C 1, deviations of Kr from K are expected to be small. Thus, considering for instance the constraint 
|1 - Kr/K\ < 0.1, one gets a(A +B) < 0.2; the limiting case K ^ Kr gives K ~ (-1 ± 1.1) x 10"^, indicating 
a perturbation with a very short range (actually, well inside the Sun, so that one should work with the interior 
Scharzschild solution). The range of allowed values for these parameters is depicted in Figures 1 and 2. 

In the limitM — » 0, Eqs. (18) yield gtt = —e^^l^ and grr = —e^^''^^. An analogy with Rosen's bimetric theory allows 
for the PPN parameter to be obtained, by interpreting this change of the metric component as due to a background 
metric rfnv [19]. Notice, however, that the vector field no longer rolls to a radial vev in the absence of a central mass, 
since this spatial symmetry is inherited from its presence, so this result should be taken with caution; this said, one 
obtains [20] a2 = (c/cg) - 1 =e2(^r-^)/'? _ 1 ^ 2{Kr-K)/^ = (A +B)a/<^, which has a radial dependency. Assuming 
a ~ 1 and considering the spin precession constraint arising from solar to ecliptic alignment measurements [17], one 
has \a2\ = |A+B|/r<4x 10"'^, implying that |A+B| <4x lO~'^rQ = 2.78 x lO^m, whererg = 6.96 x lO^misthe 
radius of the Sun. 

Further pursuing this analogy, we remark that, since there is no explicit Lorentz breaking, the speed of light remains 
equal to c. However, the speed of gravitational waves cg is shifted by an amount 



FIGURE 1. Allowed values for K and Kr [11]. 




£2 = e^(K-Kr)l£, 



1 



A+B 



(20) 



and hence it acquires a radial dependence. As stated before, this result is highly simplistic and should be taken with 
caution, since it lacks a complete treatment of gravitational radiation induced by LSB, taking into account variations 
of the bumblebee field B^u around its vev b^. 

Finally, notice that, since the radial LSB effects are exact, while the radial/temporal results are not, a direct 
comparison of these scenarios by taking the A ^ limit is not possible. 



4. AXIAL/TEMPORAL LSB 

The anisotropic LSB case is dealt with in this section. As before, we assume that the bumblebee field is stabilized at 
its vacuum expectation value, which possesses both a temporal and a spatial component; the latter is taken to lie on the 
X-axis, that is, b^i = JC^' (a,Z7,0,0). Since the radial symmetry of the Scharzschild is clearly broken, one cannot resort 
to a Birkhoff canonical Ansatz- Instead, the perturbations h^y to the flat Minkowsky metric must be obtained. To first 
order mh^iy, one has 



^00 • 

^0 = 



where time derivatives were neglected, since one assumes that v <C c. 

In order to solve the Einstein equations, one first writes the stress-energy tensor for the bumblebee field. 



Tbhv = 



^b"bPRapg^,y - b^^b^Rav - byb"Ran 



which has a vanishing trace. From the trace of the Einstein equations, one gets 

1 



R^iv — ^ 



To get the /zqo component to first order in the potential U, one writes 

^00 = - ^ [a^Rm + b^Rn + 2abRw) - 2a {aRoo + bRw) 
which, after a little algebra [11], yields the differential equation 

2 + 5a^-b^\, 

— 2_|_g^2 — ) "00,11 + «00,22 + «00,33 = . 



This admits the solution 



hoo{x,y,z) 



2M 



where ^ = (2 + 5a^)/{2 + 5a^ - b'^). 
Similarly, the ha components {i^V) obey 



(2 - b^) +4hu,u + 2hii,jj = {a^ - b^) /zoo,ii + + hooji + a^hoojj ■ 

Taking the Ansatz hu{x,y,z) = hoo{ccix, Oiy, aj,z) and, after some calculation (see Appendix I of Ref. [11]), 
obtain for the coefficients a,: 



a{ = \ 



(2-a 



2\u2 



Hence, 



with the definition 



(2 + 5a2)(2-i,2) 
hii(x,y,z) =h(M{a\x,y,z) = 



at = at 



2M 
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^Jclx^+y^ + z^ 

2{2 + 5a^-2b^)-4a^b^ 



(2 + 5a2-^2)(2_^2) • 
The hii component is now computed; a similar calculation leads to the differential equation 1 1 1] 



h 



11,22 + '^11,33 = 2M 



2^-2 n \ 2cy-y^-z^ ^2 2clx^-y^-z^ 



2 + 3/72 



'V + y2 + z2)5/2 ^^2 ^^2^2 +^2 + ^2)5/2 

indicating that the solution is a hnear combination of hm and h22- Indeed, 



hnix,y,z) = - ( 2 + 3b^ +cl]hoo{x,y,z)+2cjh22{x,y,z) ■ 
Proceeding to the off-diagonal component hio, one obtains the differential equation 

/«oi,22 + 'Joi,33 = --r-, — 9 , ,9 [hoo,22 + hoo,33+hn,22 + hn,33 — '2h22,n] 
l + a-^ + b^ 

Writing /jqi = —ab/{\ + b^){h(yQ+hn + 5hoi) leads to 

4Mcl{2clx^ - - Z^) 



5/l01,22 + 5/loi,33 = 2/j22,ll 



(c2x2+);2+z2)5/2 



and hence 5/ioi(^,)',z) = —2c2h22{x,y,z). Therefore, 



hoi = 



ab 



<^Hcl-i) . -2 

2 + 3fc2 +^^0 



(32) 



(33) 



(34) 



(35) 



l+a2-|.^2 

Finally, the /jqo component is computed to second order (see Appendix II of Ref . [11]); it can be shown that only a 



correction to the first order term /iqq appears: 
hoo = 



2cl [6 + 9fc2 + ( 1 5 + 22b^)a^] + a^b^ ( , ) 



c2(6 + 15fl2 + /,2)(2 + 3/72) 



'00 



00 



(36) 



The PPN formalism cannot be straightforwardly used to ascertain its effects, since it relies on a transformation to 
a quasi-cartesian frame of reference on which, by definition, all diagonal metric components gu are equal. However, 
some PPN-like parameters may be extracted from the results, by noticing that 



1 



1 



l-yfcco/0 



For hoo, one gets 



^{l+2k)x^+y^ + Z^ r 

\-(cl-\)cos^Q 



hoo '■ 



b^\ . 
1--12M- 



(37) 



(38) 



Since no r^2 correction appears, the PPN parameter j3 vanishes in this approach. However, ashu =/= h22 = h^i-, the same 
reasoning allows two parameters analogous to the 7 PPN parameter to be obtained: after neglecting the normalization 
with respect to hoo, one gets 



7i = H- cos 6 X 



2 + 3^2 +^0 



ab 



(l-cl)+2cl{\-cl) 



. 2. 

~ 1 + —cos^e 



i + 4)cos^e ~ 1 - ( y ) cos^e . 



(39) 



As expected, the x-axis LSB produces a stronger effect on the /in component. No clear connection can be derived 
to link 7 with ji and 72, due to the aforementioned anisotropy. However, one can take 7 to be of the same order of 
magnitude as the average of 71 and 7^, integrated over one orbit: 



1 

7- 1 - 2 (7i + 72) - 1 ^ ^ (cos^e) 



(40) 



where e is the orbit eccentricity. For low values of e, one gets 7 ~ The constraint 7 = 1 + (2. 1 ± 2.3) x 10 ^ then 
enables < 1.9x lO'^. 

A discussion concerning the anisotropy of inertia and its effect in the width of resonance lines has been presented 
as a test between Mach's principle and the Equivalence Principle [21, 22], relying on the hypothetical effect on 
the proton mass of the proximity to the galactic core. In the present scenario, we note that a radial LSB with the 



galactic core acting as source would amount to an axial LSB in a small region such as the Solar System. The bound 
Awp/wp < 3 X 10~^^, mp being the proton mass [23], can then be used to obtain 



In this contribution, the solutions of a gravity model coupled to a vector field where Lorentz symmetry is spontaneously 
broken are studied, and three different relevant scenarios were highlighted: a purely radial, temporal/radial and 
temporal/axial LSB. 

In the first case, a new black hole solution is found, exhibiting a removable singularity at a horizon of radius 
Tj = {2Mri^^y^^^^, slightly perturbed with respect to the usual Scharzschild radius r^o = 2M. This has an associated 
Hawking temperature of J = {IMtq^Y^Tq, where Jq = H/^kkeM is the usual Hawking temperature, and protects an 
intrinsic singularity at r = 0. Bounds on deviations from Kepler's law yield L < 2 x 10~^. 

The temporal/radial scenario produces a slightly perturbed metric that leads to the PPN parameters j3 w 1 — (/T + 
Kr)/M and 7 « 1 — {K + 2Kr) /M, directly proportional to the strength of the derived effect (given by K and Ky ~ K). 
Since K and are integration constants, no constraints on the physical parameters may be derived from the observed 
limits on the PPN parameters. Also, an analogy with Rosen's bimetric theory, yields the PPN parameter 7~ (A 
^ being the distance to the central body and A and B parameters ruling the temporal and radial components of the 
vector field vev. 

In the temporal/axial scenario, a breakdown of isotropy is obtained, disallowing a standard PPN analysis. However, 
the direction-dependent "PPN" parameters Ji ~ b^cos^Q /2 and 72 — cP'b^cos^Q /A may be derived, where a and b 
are respectively the temporal and x-component of the bumblebee vector vev; naturally, 71 <C T^. A crude estimative 
of the PPN parameter 7 yields 7 w b^(l — e^)/4, where e is the orbit's eccentricity. Furthermore, a comparison with 
experiments concerning the anisotropy of inertia produces the bounds < 2.4 x 10~^^ 
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